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Abstract 

We analyze the octet baryon masses and the pion/kaon-nucleon cr-terms in 
the framework of heavy baryon chiral perturbation theory. We include all terms 
up-to-and-including quadratic order in the light quark masses, m q . We develop 
a consistent scheme to estimate low-energy constants related to scalar-isoscalar 
operators in the framework of resonance exchange involving one-loop graphs. The 
pertinent low-energy constants can only be estimated up to some finite coefficients. 
Including contributions from loop graphs with intermediate spin-3/2 decuplet and 
spin-1/2 octet states and from tree graphs including scalar mesons, we use the octet 
baryon masses and the pion-nucleon cr-term to fix all but one of these coefficients. 
Physical results are insensitive to this remaining parameter. It is also demonstrated 
that two-loop corrections only modify some of the subleading low-energy constants. 
We find for the baryon mass in the chiral limit, m= 770 ±110 MeV. While the 
corrections of order rr? q are small for the nucleon, they are still sizeable for the A, 
the S and the S. Therefore a definitive statement about the convergence of three- 
flavor baryon chiral perturbation can not yet be made. The strangeness content of 
the nucleon is y = 0.21 ± 0.20. We also estimate the kaon-nucleon cr-terms and 
some two-loops contributions to the nucleon mass. 
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1 Introduction 



Despite two decades of calculations, there is still sizeable uncertainty about the quark mass 
expansion of the octet baryon masses in the framework of chiral perturbation theory, see 
e.g. 0, |, H |, |, |, 0, i |, H, |H U, HI. The penultimate goal of these calculations 



is to get additional bounds on the ratios of the light quark masses and to understand 
the success of the Gell-Mann-Okubo mass relation. Originally, it was argued that the 
leading non-analytic corrections of the type M|, where denotes the Goldstone boson 
fields 7T, K, 7], are so large that chiral perturbation theory is meaningless in this case (for 
the most extensive study see Gasser ||). In the same paper, a meson-cloud model was 
proposed which showed a faster convergence and allowed to give some bounds on the quark 
mass ratios. With the advent of a heavy baryon chiral perturbation theory as proposed by 



Jenkins and Manohar ||14| , which allows for a consistent power counting and is in spirit 
similar to the meson-cloud model, the problem of the chiral expansion of the baryon 
masses was taken up again. In particular, it was argued by Jenkins || that the inclusion 
of the decuplet in the effective field theory provides a natural framework to account for 
the large cancellations between meson loop and contact term contributions at orders 
q 2 , q 3 and (partly) q A (where q denotes a small momentum or meson mass) thus giving 
credit to the Gell-Mann-Okubo (GMO) mass relation. However, the inclusion of the 
decuplet in that approach was not done systematically and criticized in Ref. .0 Lebed 
and Luty M analyzed the baryon masses to second order in the quark masses, claiming 
that the quark mass ratios and deviations from the GMO relations are not computable, 
whereas the corrections to the Coleman-Glashow relations in the decuplet are. We have 
already pointed out in the letter M] that in that paper not all contributions at order 



m 



and m 2 In m q are accounted for. Furthermore, most of the existing papers considered 
mostly the so-called computable corrections of order m 2 q (modulo logs) or included some 
of the finite terms at this order |5], ^, || |1(J] . This, however, contradicts the spirit of chiral 
perturbation theory (CHPT) in that all terms at a given order have to be retained, see e.g. 



17], [Tj| [19| . In Ref. [|16[] we presented the results of a first calculation including all terms 



of order 0(m 2 ), where m q is a generic symbol for any one of the light quark masses m u 
In that paper, the isospin limit m u = was taken and electromagnetic corrections were 
neglected. Concerning the latter, we remark here that there is a fundamental problem 
in that the electromagnetic interactions lead to an infinite renormalization of the quark 
masses as first pointed out by Gasser and Leutwyler | 20fl . This problem has yet to be 
addressed in a systematic fashion and we thus will work in a world with a = 0, with 
a = e 2 /4iT the fine structure constant. We note that the treatment of estimating low- 
energy constants from intermediate decuplet states in Ref.jn5] was inconsistent. Here, 
we will give the correct and much more detailed treatment of the s-channel resonance 
exchanges to the low-energy constants and thus refine our previous analysis |T(| . 

Closely related to the chiral expansion of the baryon masses is the pion-nucleon a-term 
H |7|. |T[] (and also the two kaon-nucleon a-terms, a K ' N ). These a-terms measure the 

# 4 Notice that only recently a systematic effective field theory formulation for spin-3/2 fields has become 
available for the two-flavor case [tl5|l . 
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strengths of the scalar quark condensates in the proton and vanish in the chiral limit of 
zero quark masses. They are thus particularly suited to test our understanding of the 
mechanism of the spontaneous and explicit chiral symmetry breaking. Again, at present 
no clear picture concerning the quark mass expansion of these quantities has emerged, 
although it is believed that the dispersion-theoretical analysis of Gasser, Leutwyler and 
Sainio []22[ has significantly sharpened our understanding of relating a^N as extracted from 
the pion-nucleon scattering data to the expectation value of m(uu + dd) in the proton, 
where m = (m u + 7?\z)/2 denotes the average light quark mass. Of particular interest is 
the so-called strangeness content of the proton, i.e. the matrix element < p\m s ss\p >. 
It can be indirectly inferred from the analysis of the pion-nucleon cr-term. The present 
status can be summarized by the statement that this matrix element is non-vanishing 
but also not particularly large, as one deduces from the lowest order quark mass analysis 
(for a status report, see e.g. Ref.|23|]). Our aim is to show to what accuarcy one can 
at present make statements about the quark mass expansion of the a-terms. For some 
comprehensive reviews, see e.g. [[24], |25|, |26], |27j. 

The paper is organized as follows. In section we write down the effective meson- 
baryon Lagrangian, i.e. all terms which contribute to the baryon masses and the cr-terms. 
The minimal number of contact terms up-to-and-including 0(q 4 ) consists of seven terms 
of dimension two and seven terms of dimension four. These terms are accompanied by 
so-called low-energy constants (LECs), i.e. coupling constants not fixed by chiral sym- 
metry. We remark that some of the dimension four LECs can be absorbed in some of 
the dimension two ones since they amount to quark mass renormalizations of the latter. 
This ambiguity is also discussed. In section |3] we give the formulae for the fourth order 
contribution to the octet baryon masses and the pion-nucleon a-term. We also discuss 
the renormalization to render the divergent loop diagrams at order q 4 finite. This is in 
marked contrast to the calculation of the leading non-analytic corrections which are all 
finite. Section ^ deals with the estimate of the fourteen LECs from resonance exchange. 
We consider contributions from the spin-3/2 decuplet, the spin-1/2 octet of even-parity 
excited baryon resonances and from t-channel scalar meson exchange. Considering loop 
graphs with intermediate resonances leads to divergences to leading order. The corre- 
sponding renormalization constants can not be fixed from resonance parameters. This 



problem was incorrectly treated in fL6 |. However, using the baryon masses and the pion- 



nucleon er-term, we are able to fix all but one of these coefficients. We show that the 
corresponding LECs and physical observables are essentially insensitive to the choice of 
this parameter which is bounded phenomenologically. We also address the question of 
two (and higher) loop contributions to the LECs. This is mandated by the fact that 
the estimation of the LECs for the baryon masses and a-terms involves Goldstone boson 
loops (treated in a particular fashion as explained in section |J) and intermediate resonance 
propagators. Therefore, a non-trivial reodering of the chiral expansion for such graphs 
is expected. As we will show two-loop graphs only contribute to some of the subleading 
LECs at the same order as the one-loop diagrams. Most of the two-loop contributions 
can be completely absorbed in a redefinition of the one loop renormalization parameters. 
The results are given in section [| including a detailed study of the dependence on certain 
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input parameters. A short summary and outlook is given in section |]. Some lengthy 
formulae and technicalities are relegated to the appendices. 



2 Effective Lagrangian 

To perform the calculations, we make use of the effective meson-baryon Lagrangian. To 
construct the various terms, we start from the relativistic formulation which is then re- 
duced to the heavy fermion limit. This has the advantage of automatically generating all 
kinematical 1/m corrections with the correct coefficients, see e.g. appendix A of Ref. p8 



However, we do not spell out the relativistic Lagrangian but directly the heavy mass 
form emerging from it. In the extreme non-relativistic limit Jl4|, ^8[, the baryons are 
characterized by a four- velocity t> M . In this approach, there is a one-to-one correspon- 
dence between the expansion in small momenta and quark masses and the expansion in 
Goldstone boson loops, i.e. a consistent power counting scheme emerges. Our notation is 
identical to the one used in ]7| and we discuss here only the new terms necessary for the 
calculations of the masses and a-terms. The pseudoscalar Goldstone fields (0 = tt, K, rf) 
are collected in the 3x3 unimodular, unitary matrix U(x), 



t/(0)= M 2 (0) = exp{^/F 7r } 
with F n the pseudoscalar decay constant (in the chiral limit), and 

+ > - + K+ \ 
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Under SU(3) L xSU(3) i? , U{x) transforms as U U' = LUR\ with L, R E SU(3) L)jR . The 
matrix B denotes the baryon octet, 



B= S" 
V E 



J_yo i J_a 
wo 



P 
n 



\ 



(3) 



Under SU(3)l x SU(3)r, B transforms as any matter field, 

B -> B' = KBK ] , 



(4) 



with K(U, L, R) the compensator field representing an element of the conserved subgroup 
SU(3)y. The effective Lagrangian takes the form 



C 



_ Ml) r {2) A3) r (4) r {2) r (4) 
cff — *^if)B *~-(j)B *~-<j)B *~-<j)B < '~'4> < ^<j> 



(5) 



where the chiral dimension (i) counts the number of derivatives and/or meson mass inser- 
tions. Before discussing the various terms in C e g, a few words about the power counting 
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are in order. Clearly, derivatives (momenta) and meson masses count as order q. For the 
quark masses, we stick to the standard scenario, which prescribes a power q 2 to scalar and 
pseudoscalar external sources, in particular to the eigenvalues of the quark mass matrix 
M.. It would certainly be interesting to study the alternative approach in which the quark 
masses are counted as order q only. Clearly, one then would have to deal with much more 
terms at a given fixed order. We do not attempt such a 'generalized' analysis here. We 

(3) 

remark that the dimension three operators in C^g do not contribute directly to the masses 
but are needed for the renormalization of the baryon self-energy. 
The form of the lowest order meson-baryon Lagrangian is 

Cf B = iBv D B + FTriBS^u^B}) + DTi{BS^{u^ B}) , (6) 

where D ~ 3/4 and F ~ 1/2 are the two axial- vector coupling constants, = d^ + T^ is 
the chiral covariant derivative^ and denotes the covariant spin-operator a la Pauli- 
Lubanski. The meson Lagrangian + C$ is given in |19 . 



The dimension two meson-baryon Lagrangian can be written as (we only enumerate 
the terms which contribute) 

4 2 i = <B br) +E^p + 4r , ( 7) 

8=1 

(2) 

with the 0\ monomials in the fields of chiral dimension two discussed below. The explicit 
symmetry breaking terms are collected in 

£% hr) = b D Tr(B{ X+ , B}) + b F Tr(B[ X+ ,B\) + b Tr(BB) Tr( X+ ) , (8) 

i.e. it contains three low-energy constants &o,d,f- X+ — u^x 1 ^ + ux*u is proportional to 
the quark mass matrix M. = diag(m n , m^, m s ) since x = 2BAi. Here, B = ~(0\qq\0)/F 2 
is the order parameter of the spontaneous symmetry violation. We assume B ^> F n . In 
what follows, we will work to order g 4 and thus have to include terms with derivatives on 
the Goldstone boson fields. The minimal set of such terms is given by 

hOf ] = 6 1 Tr( J BK,K, J B]])+6 2 Tr(5K,K, J B}]) 

i=l,2,3,8 

+ h Tr(5K, K, B}}) + b 8 Tr(BB)Tr(u^) (9) 

with = iwdfJJw. A couple of remarks are in order. First, there are also terms of the 
type B(v ■ u) 2 B. Because of the identity (in d space-time dimensions) 



# 5 For the calculations performed here, we only need the partial derivative in 
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which holds for all the graphs to be considered later on, the respective coupling constants 
can completely be absorbed in the 61,2,3,8- Second, we heavily make use of the Cayley- 
Hamilton identity for any traceless 3x3 matrix X, 

Ti(B{X 2 , B}) + Tr{BXBX) - ^Tr(BB)Tr(X 2 ) - Tt{BX)Tt{BX) = . (11) 

This allows e.g. to eliminate the term of the type Tr(5it At )Tr(w M 5) in terms of the Ox, 
O3 and Os, i.e. 

^Tr(B[u^, K,5]]) + ^(%,{<fl}}) 

-~Tr(BB)Tr(u ll u' t ) - Tr Tr (^5) = . (12) 



The basis used in Ref.[16] was overcomplete and thus a larger number of terms was 
counted. This does, however, not directly matter for the calculations.* 6 There are also 
relativistic corrections of dimension two, i.e. 1/m corrections with fixed coefficients. These 
read 

4b C) = -^Tr(BS,[D»,{v.u,B}])-^Tr(BS»[D»,[v.u,B}}) 
2 m 2 m 

_lL Tl (BS,[v ■ u, [D», B}}) - -^TriBS^v ■ u, [Z>", B]}) 
2 m 2 m 

D 2 — 3F 2 - D 2 
H — Tr(B[v • u, [v ■ u, B]]) — Tr(BB) Ti(v ■ uv ■ u) 

24 m 12 m 

_ _L Tr(fi [Dp, [D^ B]]) + -4- MB [v ■ D, [v ■ D, B]]) 
2 m 2 m 

DF 

-Tt(B[vu,{vu,B}]) , (13) 

4 m 

where m denotes the average baryon octet mass in the chiral limit. All low-energy con- 
stants in C?b are finite. The splitting of the dimension two meson-baryon Lagrangian 
in Eq.(^) is motivated by the fact that while the first three terms appear in tree and 
loop graphs, the latter eleven only come in via loops. We remark that the LECs 6j 
{i = 0,1,2,3,8, D, F) have dimension mass -1 . Notice also that the last three terms in 
Eq.(|13|) could be absorbed in the LECs 61,2,3- For our later resonance saturation estimates, 
we prefer to keep them separately. 

There are seven independent terms contributing at dimension four, which can be de- 
duced from the over-complete basis 



i=l 



# 6 Notice that for facilitating the comparison with the results of Ref.[[l6[, we keep the numbering of the 
terms used there. 
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= d x Tr(£[ X+ , [ X+ , B}}) + d 2 Tr(B[ X+ , B}]) 

+ d 3 Tr(B{ X+ , {X+, B}}) + d A TT(B X+ )Tr(x + B) 

+ d, Tr(% + , B])Tr(x+) + 4 MHx+, B})Tr( X+ ) 

+ d 7 Ti(BB)TT(x + )Mx + )+d 8 Ti(BB)TT(x 2 + ) ■ (14) 

We choose to eliminate the c?6^term by use of the Cayley-Hamilton identity Eg . (|TTD for 
X = X + — Tr(x + )/3. The dt have dimension mass -3 . It is important to note that some of 
the di simply amount to quark mass renormalizations of some of the dimension two LECs 
||. To be specific, one can absorb the effects of d^ and d 7 in 6p and bo, respectively, as 
follows 

b F ^b F - d 5 Tr(x+) , b b - d 7 Tr(x+) . (15) 

This is a very general phenomenon of CHPT calculations in higher orders. For example, in 
7T7T scattering there are six LECs at two loop order (g 6 ) p9f , but only two new independent 
terms ~ s 3 and ~ s M^. The other four LECs make the q 4 counter terms ii (i = 1, 2, 3, 4) 
quark mass dependent. At this point, one has two options. One can either treat the 
higher order LECs as independent from the lower order ones or lump them together to 
mimimize the number of independent terms. In the latter case, one needs to refit the 
numerical values of the lower dimension LECs. We prefer to work with the first option 
and treat all the di separately from the 6j. Therefore, we have fourteen LECs not fixed by 
chiral symmetry in addition to the F and D coupling constants (subject to the constraint 
F + D = gA = 1.25) from the lowest order Lagrangian C^ B . What we have to calculate 

are all one-loop graphs with insertions from and tree graphs from We stress 

that we do not include the spin-3/2 decuplet in the effective field theory ||, but rather 
use these fields to estimate the pertinent low-energy constants (resonance saturation 
principle). We therefore strictly expand in small quark masses and external momenta 
(collectively denoted by 'q') with no recourse to large N c arguments. 

3 Baryon masses, nN a— term and renormalization 

In this section, we assemble the formulae for the quark mass squared contribution to the 
baryon masses and the a-term. We also discuss the necessary renormalization of the 
LECs di since at order q 4 the meson loop contributions to the baryon masses and the 
cx-terms are no longer finite. 

3.1 Fourth order contribution to the baryon masses 

In general, the quark mass expansion of the octet baryon masses takes the form 

m = m+Y. B q m q + Y. C i m T + E D q m 2 q + ... (16) 

q q q 

modulo logs. The coefficients B q , C q , D q are state-dependent. Furthermore, they include 
contributions proprotional to the low-energy constants enumerated in the previous sec- 
tion. Our aim is to evaluate the terms of fourth order in the chiral counting. The lower 
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order terms ~ m q and ~ m 3 / 2 for the baryon masses are standard, we use here the same 
notation as Ref. [0]. Calculating the one loop graphs shown in Fig.l and the counter terms 
~ di not shown in that figure, the q A contribution to any octet baryon mass mg can be 
written as 

,(4) _ jp ^4 , pq ^ w„ , p M P M* p Q Ml Ml Ml, 



<> = <b K + eg Ml Ml + ln(^) + ln(^f) , (17) 

X X 



with P, Q = 7r,K,i] and A the scale of dimensional regularization. Also, A x = &TrF n is 
related to the scale of chiral symmetry breaking. In fact, this is the canonical pref actor 
appearing in one loop integrals as performed here. We remark that the graphs with an 
insertion from C^ B on the baryon propagator lead to a state-dependent shift, i.e. they 
contribute to the baryon mass splittings at order q 4 . The exception to this rule is the 
term proportional to 6q- Expanded to quadratic order in the meson fields, it takes the 
form 

4 B b Tr(BB) Tr(M) -^-Bb Tr(BB) Tr(0 2 M) + £>(0 4 ) (18) 

pi 

and thus its contribution to the baryon mass can be completely absorbed in the octet 
mass in the chiral limit, 

m Ti(BB) - b Tr(BB)Tr(x+) =m Tr(BB) . (19) 

Therefore, as in the case of the g 3 calculation || [F| one needs additional information 
(like from one of the er-terms) to disentangle bo and m. The explicit form of the state- 
dependent prefactors for the nucleon, i.e. the ef_jy , is (for the LECs di the renormalized, 
finite values enter as denoted by the superscript V, see next paragraph, and we do not 
show their explicit scale dependence) 

el N = -4{4d[ + 2d\ + 4 + 3d r 8 ) , e n 2 K N = 8{4d[ - 2d r 2 -d\- 2d r 7 + 2<Q , 
ef N = 16{-d[ + cT 2 - d r 3 + d\ - d r 7 - d r 8 ) , 

4,n = 3(&d + b F + 2b - h - b 2 - b 3 - 2b 8 ) - 31 D "" '' r 



4, N = 2 ( 3 ^ - b F + 46 - 36i + b 2 - 3b 3 - 46 8 ) 



4 m 

5D 2 - 6DF + 9F 2 
6 m 



e\ N = 2b - 36! + b 2 - \h - 2b 8 - ^ 3 p* , e^ N = , e\ K N 

<J 12 m 

9 



ef* = (5D 2 -QDF + 9F 2 )(2b F + b )--(D-F) 2 (2b D + b ) 

-l(D + 3F) 2 (3b -2b D ) , 
o 

e££ = -2(5D 2 -6DF + 9F 2 )(-b D + b F ~b )-9(D-F) 2 b 

~(D + 3F) 2 (3b + 4b D ) , 

= -6u + ' e tN = ^ b D - ^b F . (20) 



The corresponding coefficients for the A, E and the H are collected in Appendix [A[ A few 
remarks are in order. First, we have e[ B = e£" B = e± B = making use of the Gell-Mann- 
Okubo relation for the pseudoscalar meson masses. In these terms, the deviations from 
the GMO limit enters only at order q 6 . Second, the e^ B could be absorbed in the ef B . We 
prefer to keep the more symmetric notation given in Eq . (^0[) . Third, we stress that the 
mixed terms ~ Mp Mq were not considered in most of the existing investigations. Notice 
that we have kept the baryon mass in the chiral limit. In the fourth order terms, it could 
simply be substituted by the physical mass (the difference being of higher order). 



3.2 Renormalization 

The fourth order contribution to the baryon masses are no longer finite. The appearance 
of these divergences is in marked contrast with the q 3 calculation which is completely 
finite (in the heavy fermion approach) . To be precise, we calculate the baryon self-energy 
Sb(w), which is related to the baryon propagator Sb(u) via 

SbH = °' „ > ( 21 ) 

p ■ v— m — Lb (a;) 

with p^ —m v fl + Z M , v ■ I <Cm and to — v ■ I. The propagator develops a pole at p = m B v, 
with ins the renormalized (physical) baryon mass, 

rriB —Ui +L B (0) . (22) 

The nucleon wave-function renormalization is determined by the residue of the propagator 
at the physical mass, 

S B (u) = , Z B = 1 + L' B (0) . (23) 

p ■ V — TUB 

To calculate the pertinent loop graphs, we use dimensional regularization and all appear- 
ing divergences can be absorbed in the LECs di 

di = d^X) + ^L. (24) 

with A the scale of dimensional regularization and 

L = ^-{ih-h lD< ^ +1 -^}- < 25) 

Here, 7# = 0.5772215 is the Euler-Mascheroni constant. The scale dependence of the 
4(A) follows from Eq.(g4|): 

<(A 2 )=<(A l) + T5 ^.n^ . (26) 
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In what follows, we set A = 1 GeV. The corresponding Tj read 
Ti = + ^b 3 + -L(7 + 15D 2 + 27 F 2 )b D 



+ 



1 —{D 2 - 3F 2 ) - | {-\{D 2 - 3F 2 )6 + 2 -DFb F ) 



72 m 



T 2 = ^ 2 + -^-(l + 24 J D 2 )6 F + ^ J DF 
4 12 96 m 

T 3 = -\b l -^b z + \(2 + 3{D 2 -3F 2 ))b D 

+ — 7 —(D 2 - 3F 2 ) - - A (--(D 2 - 3F 2 )b - 3DFb F ) 
128 m 4 V 4 / 

3, 1, 1/11 9,68 ~ A -A, 

T 4 = -bi 6 3 + - — D 2 -AF 2 ))b D 

2 18 2 V 9 4 V 27 V 

4rP 2 -3F 2 )-^ 2 -F 2 k 
8 m 4 w / 

18 2 2 l 9 3 ;F 2 ° 36m 

1, 35, 11, 7 ,35 „ 

r 7 = — h 63 b 8 -(— D 2 + F 2 ) 

4 1 108 3 18 8 128 m 27 ; 

1 / 22 34 , ,\, 1 /44 43 9 ,\, 



8V 9 3 / " 8V9 

r » - -i 6 '-i 6 '-^-i^'ir B2 + F2 »- 3B ^ 
- K-f + 5 fl2+6F2 ) 6D + 5(5-i (B2 - 3F2 »>» 

We note that for performing this renormalization, one has to include the terms 

-2i {^D 2 + F 2 ) Ti(B[v ■ D, [v ■ D, [v ■ D, B}}}) 

+i\{D 2 - 3F 2 ) Tr(5{ X+ , [v ■ D, B]}) - i^DFTr(B[ X+ , [v ■ D, B]]) 
3 13 

3 ' ->2 o r 2 



-(D 2 - 3F 2 )Tr(B{ X+ , [D», [D„B]]}) 
8 m 

+A r DFTr(B[ X+ ,[D^[D„B}}}) 
4 m 

+ l(y^ 2 + **WB[£>", [D„ B]])Tr(x + ) 



+ {|( J D 2 - 3F 2 ) - 3^D 2 b - &d(^ 2 - 3F 2 ) - QDFb F ^j j 
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+{-—DF - ±(b F ($D 2 + 9F 2 ) - lODFb D ^J } Tr(B[ X+ , [v ■ D, [v ■ D, B]]]) 
+ {-4(f^ 2 + F 2 ) — 6 (b (^D 2 + 6F 2 ) + b D (2F 2 + 2 ^D 2 ) + ADFb F ) } 



xTr(B[D-A[«-A5]])Tr(x + ) 



+ 4(10£ 2 + 18F 2 ) Tr(5[t; .£),[«• D, [D», [£>„, B]]]]) 
m 

~ C^-D 2 + 6F 2 ) Tr(B[v ■ D,[v D, [v ■ D,[v ■ D, B]]]]) (28) 

777, \ o / 

in the effective Lagrangian (an overall factor LF~ 2 has been scaled out and the finite 
pieces related to these operators can be neglected). These, however, do not directly 
contribute to the masses and cr-terms and are therefore not listed in section 3.1. It is 
obvious that the dimension two terms also enter the renormalization procedure and thus 
the finite constants bi appear in the IY This is a general feature of any renormalization 
beyond 0(q 3 ) in heavy baryon CHPT. 

3.3 Pion/kaon— nucleon a— terms and strangeness content 

Further information on the scalar sector of baryon CHPT is given by the scalar form 
factors or cr-terms which measure the strength of the various matrix elements m q qq (q = 
u, d, s) in the proton. In a mass-independent renormalization scheme, one can define the 
following renormalization-group invariant quantities: 

c-kN (t) = m < p \uu + dd\p > , 

°"SvW = ^(m + m s ) <p'\uu + ss\p> , 

a$ N (t) = ^(m + m s ) < p' | - uu + 2dd + ss\p > , (29) 

with t — (p' — p) 2 the invariant momentum transfer squared. The explicit form of the t- 
dependent cr-terms is rather lengthy and not very instructive. Here, we just discuss some 
general aspects and refer the reader to Ref. |3(J for details. The most striking new feature 



compared to the 0(q 3 ) analysis |7|, |28], |36| is the appearance of t-dependent divergences. 
To be more specific, consider a typical diagram as shown in Fig. 2. It is most convenient 



to calculate these diagrams in the Breit-frame p8 j. For the case at hand, we get 



M 2 a --- -(AM 2 - t) In -f - (2M 2 - t)* J arcsin 



16tt 2 V a 2 2 V a ' A 2 v a ' y/i V 2M C 

where we have suppressed the pertinent Clebsch-Gordan coefficient and M a stands for a 
Goldstone boson mass in the SU(3) flavor basis. The renormalization of this diagram and 
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the other contributing to the cx-terms is somewhat tricky and is described in appendix [B|. 
Note that when combining all the terms from the various Feynman graphs (after renormal- 
ization), one is indeed left with terms proportional to const-t. In fact, these t-dependent 
coefficients can not easily be estimated. Also, they are different for a W N{t) and cr^^(t), 
respectively. The empirical information on a^{2Ml) — 0^(0) = 15MeV is thus not 
sufficient to make predictions for the two KN cr-terms shifts, a^'^ (2M] ( ) — cr^\o). We 
will therefore only give some rough estimates by simply parametrizing these t-dependent 
pieces in terms of coefficients which we can estimate to be of order one. Note that to order 
q 3 , the shifts from t — to the respective Chang-D ashen points are finite and free of any 
LEC. In view of this, the nice result found in Ref.|7| for Act^n — 15MeV by including 
the A(1232) explicitely as intermediate state must be considered accidental. 

At t = 0, the corresponding formulae simplify considerably. The fourth order contri- 
bution to ova^O) takes the form 



4>) = Ml 



Ml . Ml. P n Ml . ,Mf 



A 2 v A 2 6 ' a A 2 v A 2 

X X 



(31) 



with the ef^ are related to the efj^, Eq.(|(]), and they are collected in Appendix 0. 
These relations are, of course, a consequence of the Feynman-Hellmann theorem which 
states that 

^(0)=m^ . (32) 
am 

The derivatives with respect to the quark masses can be converted into derivatives with 
respect to the Goldstone boson masses squared, for example 

Eq.(^) can therefore be evaluated easily. By a similar reasoning, one can deduce the 
strangeness content of the proton, i.e the strength of the matrix element < p\m s ss\p >, 

i _ , dm N , . 

< p\m s ss\p >= m s — , (34) 



dm* 



and the strangeness fraction y, 

2 <p\ss\P> Ml ( 2 1 -A" 1 dm N 

V= ; = — Mz\ m,- . 35 

y <p\uu + dd\p> ct tN (0)\ k 2 V s dm s V ; 

The fourth order contribution to y, denoted by can be decomposed in the standard 
fashion (for the lower orders, see |7|) 



p »,2 p M P , , m p. pq m p , ,Mq 
< y M 2 P + el y -jf ln(^f ) + e P , Q y ln(^ 

X X 



(36) 



and the are again related to the , see appendix £]. The more lengthy expressions 
for the two KN a-terms can be found in Ref . |30H . This completes the formalism to study 
the scalar sector of baryon CHPT to fourth order in the meson masses. 
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4 Low— energy constants from resonance saturation 



Clearly, we are not able to fix all the low-energy constants appearing in from 
data, even if we would resort to large N c arguments. We will therefore use the principle 
of resonance saturation to estimate these constants. This works very accurately in the 
meson sector |31], [32], |33|. In the baryon case, one has to account for excitations of meson 
(R) and baryon (N*) resonances. One writes down the effective Lagrangian with these 
resonances chirally coupled to the Goldstone bosons and the baryon octet, calculates 
the Feynman diagrams pertinent to the process under consideration and, finally, lets the 
resonance masses go to infinity (with fixed ratios of coupling constants to masses). This 
generates higher order terms in the effective meson-baryon Lagrangian with coefficients 
expressed in terms of a few known resonance parameters. Symbolically, we can write 

C eS [U,B,R,N*]^£, eS [U,B] . (37) 

It is important to stress that only after integrating out the heavy degrees of freedom from 
the effective field theory, one is allowed to perform the heavy mass limit for the ground- 
state baryon octet. Here, there are two relevant types of contributions. One comes 
from the excitation of intermediate baryon resonances, in particular from the spin-3/2 
decuplet states. Concerning the higher baryon resonances, only the parity-even spin- 
1/2 octet including the Roper iV*(1440) plays some role. In addition, there is t-channel 



scalar and vector meson exchange. It was already shown in Ref. [16[] that the vectors do 
not come in at order q A . It is important to stress that for the resonance contribution to 
the baryon masses, one has to involve Goldstone boson loops. This is different from the 
normal situation like e.g. in form factors or scattering processes. One could argue that 
scalar meson exchange alone should provide the necessary strength of the scalar-isoscalar 
operators ~Tr(x + ). However, taking the phenomenological bounds on the scalar masses 
and coupling constants, these contributions are too small to explain e.g. the strengths of 
the symmetry breakers &o,d,f which have been determined to order q 3 in 0. To be specific, 
if one performs the calculation in SU(2) and demands that the LEC c\, which can be fixed 
from the ttN cr-term, to be saturated by scalar meson exchange only, one would need a 
mass to coupling constant ratio of M s /g s = 220 MeV. Here, g s refers to the SU(2) coupling 
g^N^N S. Such a number can not be obtained from standard scalar meson masses and 
their couplings to the nucleon (baryons). In Ref.[^] it was, however, shown, that one-loop 
graph with intermediate decuplet states can effectively produce such couplings. However, 
treating the resonances relativistically leads to four major complications. 

o First, terms arise which are non-analytic in the meson masses. Clearly, to avoid 
any double counting and to be consistent with the requirements from analytic- 
ity, one should only consider the analytic terms in the meson masses generated 
by such loop diagrams. In fact, standard tree level resonance excitation generates 
the dimension two operators fei,...^ in the conventional fashion. Calculating tad- 
pole diagrams to order q 4 with these hi generated from resonance exchange exactly 
produces the non-analytic terms ~ M? In M? which also arise from the one-loop 
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graphs with intermediate decuplet states. This is why we have to dismiss such 
non-analytic contributions in the Goldstone boson masses in our estimation of the 
scalar-isoscalar LECs from loop graphs. We remark that this kind of equivalence 
has been already discussed for the SU(2) LEC C2 which enter the calculation of the 
nucleon electromagnetic polarizabilities to order q 4 , see [|34| [|35|]. 

o Second, to the order we are working, even the remaining analytic pieces are divergent 



(this was incorrectly treated in ref. [[UJ). Therefore, we can only determine the 
analytic resonance contribution up to three renormalization constants, two related 
to the decuplet (/5a, 5a) and one to the Roper-octet (Pr). We solve this problem 
in the following manner. In addition to these three constants, we have as further 
unknowns the scalar couplings F$ and D$ and the octet baryon mass in the chiral 
limit. Fitting the four octet baryon masses and 0^(0) allows us to express all 
LECs and physical observables in terms of (3r solely. Fortunately, all quantities 
are very insensitive to the choice of this parameter which we can bound by some 
phenomeno logical argument. 

o Third, since the baryon excitations are treated relativistically as explained above, 
there is no more strict power counting |[6] and thus we must address the question 



of two-loop contributions from diagrams with baryon excitations. To leading order 
in the resonance mass expansion, they give rise to the same divergences as the one 
loop graphs and the corresponding corrections can be completely absorbed in a 
redefinition of the parameters appearing at one loop. We also show that there is 
an additional divergence at order in the Roper-octet graphs which modifies the 
finite one-loop parameter Sr. However, the contributions to the baryon masses and 
the cr-terms ~ Sr are very small and thus one does not need to know this parameter 
accurately. Using dimensional analysis, the two-loop corrections to (3r are found to 
be modest. For details, see paragraph [L5| 



o Fourth, the renormalization of the cx-terms is somewhat complicated. The di- 
vergences of the scalar form factors are renormalized through terms of the form 
dC/dm q . In general, one has dC/dm u 7^ dC/dm d 7^ dC/dm s . If one calculates in 
the isospin limit m u = one can not disentangle the derivatives with respect to 
m u and rrid any more. For the general renormalization procedure, one has to work 
with m u 7^ rrid- Applying resonance saturation, one effectively determines the LECs 
of the Lagrangian but not the ones related to the operators uu, dd and ss. By use of 
the Feynman-Hellmann theorem, the LECs related to the operators qq can be fixed 
from the ones which appear in the calculation of the baryon masses. This general 



theorem is proven in ref.|30[. It means that here we are only able to give the full 
order q A contribution to the pion-nucleon cr-term at t = and can only determine 
Urn (0) up some renormalization constant (since we work in the isospin limit). We 
also point out that we have not yet been able to generalize these methods to the 
momentum-dependent parts of the scalar form factors since there integrals appear 
which we can not work out analytically (for details, see ref.pO|). 
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Let us now consider the one-loop contributions from the decuplet, the | + -octet and 
from the scalar mesons, in order. 

4.1 Decuplet contribution to the low— energy constants 

Consider first the decuplet contribution. We treat these field relativistically and only 
at the last stage let the mass become very large. The pertinent interaction Lagrangian 
between the spin-3/2 fields (denoted by A), the baryon octet and the Goldstone bosons 
reads 

£ab* = \ { A^ c e^(z) K)i Bl e cl] + B) ( u v )t e, M (z) A£ 6c } , (38) 

with 'a,b,...,j' are SU(3)/ indices and the coupling constant C = 1.2... 1.8 can be 
determined from the decays A — > Bit. The Dirac matrix operator G Ml/ (Z) is given by 

Q^(Z)=g^-(z+l)w v . (39) 



2, 

For the off-shell parameter Z, we use Z = —0.3 from the determination of the A contri- 



bution to the Ti N scattering volume 033 |37| . That is also consistent with recent studies of 



A (1232) contributions to the nucleons' electromagnetic polarizabilities |38] and to thresh- 



old pion photo- and electroproduction |3"5fl . For the processes to be discussed, we only 
need the expanded form of u^, 



lu. 



The propagator of the spin-3 / 2 fields has the form 



a = -TT« + O(0 2 ) . (40) 



° M = ~ l p^ni - 3 W " 3< + 3m A > ' (41) 



with m& = 1.38 GeV the average decuplet mass. We have to evaluate diagrams as shown 
in Fig. 3. With the labeling of the momenta as in that figure, this leads to 

t m zl^l f d ' k (k-py&« P (z)G^(k)e, u (z)(k- P y 

A[P) 2F2 J (2tt) 4 {k-pf-Ml + ie ' 1 ' 

where the pertinent Clebsch-Gordan coefficient has been omitted and M a as defined after 
Eq.([H]). This integral is evaluated on the mass-shell of the external baryons, i.e at jj =m, 
and split into various contributions according to the power of momenta in the numerator 
and the number of propagators. Each such term is then is expanded in powers of Goldstone 
boson masses up-to-and-including O(M^) (modulo logs). Only then the large mass limit 
of the decuplet is taken. This then gives the contribution to the various LECs bi and di 
modulo some renormalization constants. To be specific, consider the scalar integral to 
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which all contributions arising in evaluating I a can be reduced 



d 4 k 



(2vr) 4 [{k - p) 2 - Ml + ie][k 2 -m\ + is] 



2L + 



16tt 2 



9i 9 ° 2 

m A i / m A\ 1/2 ° 2 \ i / m A~ m - 



m 



A 2 



72K- m ) ln (" 



+ 



+ 



2 / 2 

m {m A — m ) 



m 

l A\ , / 2 



A 2 



2 „ ,m 2 s , o o 2 x . ,m\— m 



m —m A ln(—— ) + (m A + m ) ln(- 



° 2 ^V A 2 y ' V '" A ' "" y '" v A 2 



Mi 



2(m A — m 



m A — m 



-m -2m A ln(— f ) +4m A ln( A 
A 2 

m 2 ..„. M 2 



r> o 2 

mi — m . 



A 2 



(m 2 A — m 



(43) 



In total, the M 2 \nM 2 terms cancel in I A . As explained before, the remaining non- 
analytic pieces ~ M 4 In M 2 appearing in these integrals should not be retained, see also 
below. We can now work out the complete integral Ia(J> =wi) and find at the scale A = m A 



j analytic 



-1C 2 
6F 2 



-{2Z + l)m A M 2 a - (AZ 2 + 5Z 



I m A 



Mt 2L + 



(44) 



where the ellipsis stands for constant and divergent pieces which have to be taken care of 
by standard mass renormalization terms of the type 5m A Tr(AA), see e.g. Refs.jT], |36 



and for sub leading terms in the 1 / m A -expansion at orders M 2 and M 4 . If one were to use 
the subleading finite pieces to estimate the LECs, the resulting values would be much too 
small. One notices that in this relativistic treatment, the dimension two and four LECs 
are no longer finite (as in the heavy baryon approach) [|36[j . The structure of Eq.(f44|) 
shows that the first and second term in the curly brackets will contribute to the 6o,d,f, 
and the dj, respectively, after renormalization 



ranalytic 



-1C 2 1 



2F 2 16tt 2 



-m A M 2 a a ( A ] + 



1 



m A 



M„ 4 c4 4) 



(45) 



where af' 4 ' are the finite renormalization constants related to the divergences appearing 
at order M 2 and M 4 , respectively. The 61,2,3,8 are calculated from standard pion-nucleon 
scattering tree graphs with intermediate decuplet states (for details, see ref.JJlJ). As an 
important check, feeding these back into the tadpole diagrams one finds non-analytic 
pieces which exactly agree with the ones dropped from the integral I a- These have the 
generic form 

M 4 M 2 

(2Z 2 + Z-1 )^L\n4r ■ (46) 

A 2 



7-non— analytic 
1 A 
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This proves the consistency of our procedure. Evaluating now the pertinent Clebsch- 
Gordan coefficients and defining 



e A = -^— (2Z 2 + Z-1) , p A = -f^- 2 mAa g> , S A = - C \ — ag } , (47) 
12mA 128ir z l i ^ 2bbTT z r^ m A 



we arrive at the decuplet contribution to the various LECs 

^ = o ^ A ' = _ /^A , bp = - Pa , 
o 6 

lA ^ i A lA ^ , A 3 

Oi = — r_- e A , o 2 - e A , o 3 - -4 e A , o 8 - 2 eA ' 

]_ A 1 A 1 A ^ 

^1 = - #A , d 2 = ^a , d 3 = —- 5 A , d 4 = -8a , 
df = — 6 A , d$ = — 6a , d£ = - 6 A , (4* 
where Pa and 5a have to be determined as described below. 



4.2 l/2 + — octet contribution to the low— energy constants 



The next multiplet of excited states is the octet of even-parity Roper-like spin-1/2 fields. 
While it was argued in Ref . || that these play no role, a more recent study seems to indicate 
that one can not completely neglect contributions from these states to e.g. the decuplet 
magnetic moments |40| . It is thus important to investigate the possible contribution of 
these baryon resonances to the LECs. The octet consists of the A/"* (1440), the E*(1660), 
the A* (1600) and the S*. For the mass of the latter, we use the Gell-Mann-Okubo formula, 



- m\* 
2 



- my,* 
2 



m N , 



1.79 GeV 



so that we find for the average Roper octet mass 

m R = 1.63 GeV 



(49) 



(50) 



which is 490 MeV above the mean of the ground state octet. Notice that we denote the 
spin-l/2 + octet by 'i?'.^ The effective Lagrangian of the Roper octet coupled to the 
ground state baryons takes the form 



C e g(B, R) 

C (R) 

Ant(-B, R) 



£>o(R) + £int(-B, R) , 
iTr(Ry»[D^R])-m R Tr(RR) , 
D R 



+ 



4 
Fr 
4 



Tr(Rj^ 5 { Ufl , B}) + h.c 
Tr(R^ l5 [u^B))+h.c. 



(51) 



# 7 This should not be confused with the same symbol denoting any meson resonance. 
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The numerical values of the axial-vector coupling constants D R and F R are determined in 
appendix [D| (generalizing the results for the two-flavor case 0]). With these definitions, 
we evaluate the same type of graph as shown in Fig. 3 (dropping again the Clebsch-Gordan 
coefficients), 



l! l 1 (2vr) 4 [(p - A;) 2 - M 2 + ie][£; 2 - m| + ie] ' 



which by inspection leads to the same kind of terms as in Eq . (f42|) , i.e. with k n (n = 
0,1,2,3) in the numerator combined with one or two propagators. Expanding these 
terms again first in powers of the Goldstone boson mass M a and then in powers of l/m R , 
one finds at the scale A = m R that the terms of order M 2 diverge whereas the ones of 
order M 4 are finite. After appropriate renormalization, we have 

=m) = {m R Ml a% + ^ M a 4 } . (53) 

As before, the contributions to the &i, 2,3,8 are calculated from tree-level graphs with 



intermediate Roper states (for details, see ref.|Hj). Evaluating now the pertinent Clebsch- 
Gordan coefficients and defining 

e R - ~Q > Pr = ocr 2 M m R a< R > $R = mo/i 2P2 ~ ' ( 54 ) 

Sm R 2bbn l r£ 10247t z J h ^ m R 

we arrive at spin-l/2 + octet contribution to the various LECs 

b$ = (fD R + F*)(3 R , b R = ±(3F*-D R )(3 R , b R = ~D R F R p R , 
bf = -(D 2 R -3F 2 R ) e f, b R = D R F R e R , b R = , b R = 1 -D 2 R e Rl 



< = -M-3^)|, d* = -±D R F R 8 R , d R =-(D R -3F*) S -f 



d R =(Dl-3F^ 6 f,df = D R F R ^,d R = (f 6 Dl + F^ U f 6 , 
1 47 



d§ = l ( Y Dr + Fr) 5r ' (55) 
in terms of the parameter f3 R to be fixed below. 

4.3 Scalar meson contribution to the low— energy constants 



Meson-exchange in the t-channel can also contribute to the LECs as pointed out in Jl6| 
Denoting by S and S\ the scalar octet and singlet with M s ~ M Sl ~ 1 GeV, respectively, 
the lowest order effective Lagrangian coupling the scalars to the Goldstone bosons and to 
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the baryon octet reads 



C cS (U, B, S) = C (S) + C int (U, S) + C int (B, S) , 

C (S) = ^Tr (d^S - M 2 S S 2 ) + X - Tr(d„S 1 d"S 1 - Mg t S 2 ) , 

C int (U, S) = c d Tr(Su, J ,u fl ) + c m Tr(Sx+) + c d Si Tr(w M w M ) + c m Si Tr(x+) , 
C mt (B, S) = D s Ti(B{S, B}) + F s Tr(B[S, B]) + D Sl Si Tr(BB) (56) 

where the coupling constants Ds, Fs and Ds 1 have chiral dimension zero. For the cou- 



plings of the scalars to the Goldstone bosons, we use the notation of |yj an d the pa- 
rameters determined therein, i.e. \c d \ = 42 MeV, \c m \ = 32 MeV with c d c m > 0. The 
singlet couplings can be related to these in the large-iV c limit via c d , m = ±Cd,m/ V3 with 
c d c m > 0. Similarly, we have Dg 1 = ±2D s /\/3. In fact, it is very difficult to pin down the 
couplings Ds, Fs and /^.p] We will therefore leave Fg and Dg as free parameters and 



determine them as described below. 

Evaluating now the diagrams shown in Fig.4, which are scattering graphs and direct 
couplings of the scalar mesons to the operator x+, an d using the lowest order effective 
Lagrangian Eq. ([5*B]) , one finds only contributions to Cl B - Contributions to the LECs d{ 
arise if one accounts for higher dimension operators in £ int (-B, 5"). This would amount 
to a proliferation of unknown coupling constants. We have refrained from including such 
terms based on the observation that the scalar contributon to the dimension two LECs is 
rather small [ 16| . This amounts to 



with 



s 2 

bf = ^C d D s ,b s 2 =C d F s , 6f = l -C d D s , b§ = -hj d D s + C d D Sl , 



<^m,d = ' ^m,d = -TJT ' V* 



Mg ' M Sl 



We have now assembled all pieces to discuss the numerical values of the various LECs. 



4.4 Determination of the low— energy constants 

For pinning down the numerical values of the various LECs, we need to know the values 
of the parameters Pa, 8a, Pr, Fg, F) s and of the baryon octet mass in the chiral limit, 
m. We will later show that it is mandatory to keep all these parameters, i.e. using only 
the decuplet and/or the scalar mesons to estimate the LECs does not lead to a consistent 
description of the baryon masses and a-terms. We determine these parameters as follows. 

# 8 In what follows, we neglect the singlet fields. There effects can be absorbed in the pertinent coupling 
constant of the octet fields. 
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Within the framework of resonance exchange saturation of the LECs , the baryon masses 
take the form 



m B =m+mf + X B + (3 D B + 6 D S B + eD e B + D S B , (59) 

m B - 

where we have lumped the 0(m q ) and the 0(m^) corrections togther (in the constants 
D B and D B , respectively). We remark that whenever possible, we have substituted the 
octet mass in the chiral limit by the corresponding physical mass since these differences 
are of higher order. The A contributions are the \ jm insertions from C^ B . Similarly, the 
(3 and e terms stem from tadpole graphs with insertions proportional to the low-energy 
constants &o,d,f and b i: respectively. The 5 terms subsume the contributions from C^g, 
these are proportional to the low-energy constants di. In these last three terms, we have 
abbreviated 

(3D B = P A D BA + f3 R D BtR , W B = S A D BA + 6 R D 5 BtR , 
eDl = e A D BA + e R D B>R . (60) 

Finally, the terms of the type D B are the scalar meson contributions to the mass (which 
could also be absorbed in the coefficients (3 and e, cf. Eq.(^)). The numerical values 
of the Xb, D b ,D 5 b , and D e B are given in table 1, using F = 0.5, D = 0.75, D R = 0.60, 
F R = 0.11 and F n = 100 MeV. We see that the dominant terms at Oimi) are indeed 

the tadpole graphs with an insertion from C^ B ^ (this holds for the masses but not for 
cr^jv(O)). These numbers are different form the ones given in Ref.[|IB| due to the different 



definitions of the coefficients /3, S and e and because of the Roper-octet contribution. The 
D S B are (in GeV) 

D S N = -0.013 D 5 + 0.057 F s , D'l = —0.042 D s — 0.019 F s , 

£>f = +0.052 D s + 0.019 F s , D§ = -0.045 D s - 0.076 F s , (61) 



where we have used the results of [^TJ for the scalar couplings and Ms = 0.983 GeV. The 
four equations for the octet masses 

m B = m B [m,/3 A ,8 A , (3 Rl F Sl D s \ , (62) 

allow to pin down four parameters, we choose 5a , (3 A , Fs and Ds- A similar equation can 
be worked out for the pion-nucleon a-term, 

a^(0) = a^(0) + —X a + (3 A D^+/3 R D^ + 5 A D 5 ^ + 5 R D 5 f 

m N 

+e A + e R D%* + D s D° s + F s , (63) 
or numerically with cr n ^(0) = 45 MeV 

0.077 = -0.244 (3 A + 0.020 S A - 0.056 (3 R + 0.001 D s - 0.002F 5 , (64) 
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in appropriate units of GeV. We end up with 5 equations for six unknown parameters. 
Choosing (3 R as the undetermined one, this leads to 



rh = +0.7673 + 0.3737/fe , (3 A 
5 A = -0.1387- 0.0232 , D s 
F s = -0.8974- 12.715 /3 R . 



-0.3321 - 0.1 126 f3 R 
■■ -5.8244 + 1.4484 (3 R 



(65) 



Units are GeV for m, GeV 1 for /3a,Pr and GeV 3 for 5a whereas F s and D s are 
dimensionless. The prefactors have the appropriate dimensions, for example m [GeV] = 
0.767 [GeV] +0.374 [GeV 2 ]/?/? [GeV -1 ]. We can now give a bound on (3 R based on the 
assumption that the ratio of the decuplet contribution and of the Roper octet one is 
of the same order for all LECs. Since the contributions to the 6j do not contain any 
unknown parameters, we can use their ratio to estimate the corresponding one from the 
Roper-octet and we thus conclude 



d b,r @R 



D b _aPa 



D B,R e R 



from which we derive the (soft) bound 

\P R \ < 0.1 GeV -1 . 



(66) 



(67) 



Consequently, the scalar couplings Fs and Ds are falling into the ranges — 2 < Fs < 0.5 
and —6.0 < Ds < —5.7. It would be interesting to have some phenomenological bounds 
on these couplings as a check of the consistency of our procedure. 
Finally, we can now express the LECs in terms of the parameter /3 R , 

-0.316 - 0.536 p R , 



bo 


= -0.606 -0.227/3 R 


, b D = 


+0.079 + 0.014 (5 R 


, b F 


bi 


= -0.004 + 0.024 (3 R 


,b 2 = 


-0.187- 0.421 p R 




h 


= +0.018 + 0.024 (3 R 


,b 8 = 


-0.109 - 0.032 p R 




di 


= +0.008 + 0.001 (3 R 


, d 2 = 


+0.035 + 0.006 p R 




d 3 


= +0.069 + 0.012 (3 R 


, di = 


-0.077 - 0.013 p R 




d 5 


= -0.050 - 0.008 


, d 7 = 


-0.018 - 0.003 


, ^8 



-0.103- 0.017 p R , (68) 

with the bi in GeV -1 and the di in GeV -3 and the coefficients have corresponding di- 
mensions. In table 2, we compare the LECs of the dimension two meson-baryon effective 
Lagrangian with the values previously determined from KN scattering data. We have 
transformed the results of Ref. 



into our notation. As can be seen from table 2, most 
(but not all) coefficients agree in sign and magnitude. We remark that the procedure 
used in involves the summation of arbitrary high orders via a Lippmann-Schwinger 
equation and is thus afflicted with some uncertainty not controled within CHPT.^| For 
comparison, we also give the values of 6o,d,f from the q 3 calculation [f/J. 

# 9 These values have recently been refitted taking account also rj and kaon photoproduction data p3| . 
The resulting values are somewhat different from the ones given in ref. |42| . 
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We end this paragraph with some comments on estimating the LECs from the decuplet 
or decuplet and Roper-octet alone. Consider first the case of the decuplet only. We then 
have three parameters determined by a least-square fit to the four masses. One finds 
m= 620 MeV, p A = -0.551 GeV" 1 and 5 A = 0.975 GeV- 3 and the deviations of the fitted 
from the physical masses range from -154 to 121 MeV. This means that in particular the 
deviation from the Gell-Mann-Okubo relation, 

A G mo = |(3m A + m s - 2m N - 2m s ) , (69) 

which empirically is about 6.5 MeV, is very large (of the order of 200 MeV). The LECs 
take values which are considerably larger than the ones given above (cf. table 2 and 
Eq.(^)). Furthermore, the a-term is completely fixed and turns out to be 0^(0) = 89 
MeV, considerably larger than the empirical value. One could try to remedy the situation 
by adding the Roper-octet, thus having one more free parameter at ones disposal. One 
can then solve the linear system of equations for the four masses in terms of the four 
parameters. This does, however, not lead to sensible results. One finds [3r ~ 10 GeV -1 , 
much larger than the bound given in Eq . (|67D . This leads to absurd results like m~ 
3.9 GeV. It is therefore mandatory to include the scalar meson exchange to consistently 
describe the scalar sector of three flavor baryon CHPT. 



4.5 Two— loop contributions to the LECs 

As shown, the estimation of the LECs entering the baryon masses involves Goldstone 
boson loop diagrams, compare Fig. 3. With such graphs, one encounters a new mass scale 
which is non-vanishing in the chiral limit, m R — as — > 0. Therefore, a strict 
one-to-one correspondence between the expansion in small momenta and the number of 
Goldstone boson loops is no longer guaranteed, as it is known from the study of relativistic 
baryon CHPT [[36]. A similar situation arises for the calculation of the deviations from 



Dashen's theorem, which involves photons coupled to heavy (axial-)vector mesons in the 



loops, see e.g. Refs.|44]]. Again, the new mass scale (here the vector meson mass in the 
chiral limit) spoils the strict power counting. Indeed, it has been shown that two-loop 
graphs modify the leading order results to the deviations form Dashen's theorem ||45|| . We 
therefore also have to address this issue in our context. 

The possible two-loop graphs to be considered are shown in Fig. 5. Here, we con- 
centrate on the diagrams with the maximal number of resonance propagators and the 
simplest BB*cf) coupling via the chiral connection r fl . This means for example in diagram 
(e) one of the intermediate propagators is a nucleon and the other one refers to an excited 
intermediate spin-3/2 or spin-1/2 state. For all these graphs, all possible combinations of 
decuplet, Roper-octet and nucleon propagators should be considered. From these, we only 
take the leading ones in the large mass limit subject to the constraint that the vertices 
involved are constructed from the chiral connection. The corresponding contributions 
from the axial-vector couplings ~ involve other coupling constants and are therefore 
algebraically independent from the ones considered here. In fact, to simplify the algebra, 
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we will focus on the case where all resonance propagators refer to the Roper-octet. A 
more detailed account of these calculations is given in J3U . 

Consider first the tadpole-type diagrams like in Fig. 5a. Without Clebsch-Gordan 
coefficients and other prefactors, it takes the form 
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dH 



d A k 
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32tt 2 m\ 



(70) 



where we have only retained the leading terms in the expansion in the resonance mass 
rriR as indicated by the ellipsis, i.e. we have approximated m R — m~ rriR and so on. We 
remark that the disturbing term ~ M 2 In M 2 cancels in the sum of graphs 5a, its partner 
with the tadpole at the other vertex and graph 5f. Taking into account the pertinent 
counter term graphs to renormalize the infinities connected to the tadpole, we find 
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+ O(Mt) 



(71) 



The first term in Eq. (|7T|) can be completely absorbed in the renormalization constant 
at R and the second one leads to a new renormalization constant ot R . Its finite piece 
modifies the value of the finite dimension four one-loop contribution Sr. Simarly, for 
graphs with intermediate decuplet states, to leading order in the resonace mass on can 
absorb the two-loop contribution entirely in a redefinition of the constants ot^ and at£ . 
Furthermore, if one works out the first finite contribution at orders M 2,4 , one finds that 
these are numerically much smaller than the corresponding pieces from the one loop 
diagram. Similarly, we work out the contribution from the graphs of the type 5d. The 
integral stripped off all prefactors and expanded in powers of the resonance mass takes 
the form 
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(72) 



Again, standard renormalization has to be performed and one is left with a contribution 



to at £ and one to a^' . The latter can be made to vanish if one sets A = uir. 

The other diagrams shown in Fig. 5 can not be given in closed analytical form. To get 
an estimate about their contributions, we perform asymptotic expansions in the external 
momentum p making use of the formalism developed in Ref.|46[ (and references therein). 
For external momenta (here p =m) below the first threshold (here p t h r =m +2M n ), one 
can expand around p = to leading order. It is straightforward to show PDf that this 



„(4) 
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procedure is sufficient to estimate the leading order terms in the large resonance mass 
expansion (here m R ). Specifically, we have 



in d space-time dimensions. The explicit formulae for the various contributions are very 
lengthy and can be found in Here, we just show a typical result after picking out the 
leading terms in the m R expansion, 



< 5c 



d m% + e 2 m% [M 2 a + Ml) + e 3 m R (M* + M fe 4 ) + e 4 m R M 2 a Ml + ... , (74) 



with e = d — 4 and the coefficients e« contain divergences. After removing the divergent 
pieces ~ 1/e and 1/e 2 , we are essentially left with contributions to the renormalization 
constants a^' 4 ^ and a R ' 4 \ i.e. the two-loop effects can be completely absorbed in pertinent 
redefinitions. Only the finite one-loop constant 5 R is modified. The corresponding tree 
graph contributions to the baryon masses and a-terms at order from C^ B are, however, 
very small and thus an accurate knowledge of this coupling is not needed. From this we 
conjecture that a similar mechanism is also operative at higher orders and that our one- 
loop approach to estimate the baryon resonance excitations to the various LECs is a 
consistent procedure. 



4.6 Some two— loop contributions to the baryon masses 

Having worked out the two-loop contributions to the LECs 6, and di in paragraph we 
can use the same machinery to estimate some typical two-loop, i.e. order g 5 , contributions 
to the baryon masses. This should only considered indicative since we do not attempt 
a full C(g 5 ) calculation here, which besides all two-loop diagrams would also involve 
one-loop graphs with exactly one insertion from C^g. 

Consider first the tadpole-type graph like in Fig.5a (and its partner with the tadpole 
on the other side) but with the essential difference that the intermediate propagator refers 
to a groundstate spin-1/2 state in the heavy baryon formalism. Using the appropriate 
Feynman rules, the momentum-space integral It is 

d ' q f d ' k ' ' ' -(S ■ k) 2 . (75) 



(2tt) 4 J (2tt) 4 q 2 - Ml k 2 - M 2 v (p - k) 

To leading order, we can neglect the baryon off-shell momentum. I t then takes the simple 
form 



^ = 7 f(Ma)M 2 [ ,f* * , =^M 2 a M { 



4 J y '" UJ J (2tt) 4 [k 2 - M 2 ] v ■ k 32tt "~ a " b 



1 , M 2 

2i + i^ ln # 



(76) 



As expected from the power counting, this contributions starts at order M|. Renor- 
malizing the divergence and restoring the appropriate prefactors, we have the following 
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contribution to the nucleon mass from these tadpole-type graphs 

45F 2 + 5D 2 - 6DF - , Ml 



( 5, a) _ ^ ~ U ^ J m 5 In ^ (77) 

where we have set for simplification M n = and = M v . Similarly, we can work out 
the contribution from the analog of fig. 5f with heavy baryon propagators, 

f5 f) 45F 2 + 17D 2 - 30DF - Mt 

™>n Ml\n—f . 78 

N 768tt 3 F4 x A 2 V ; 

The "double-hump" graph (compare fig. 5d) vanishes to this order, 77% = 0- Eqs. (|77| , |78D 
will be used in the next section to get an order of magnitude estimate of the two-loop 
corrections to the nucleon mass. 



5 Results and discussion 

5.1 Results for the central values of resonance parameters 

In this section, we discuss the results for our central values of parameters, setting f3 R = 
(varying (3r within its allowed range, —0.1 < (3r < 0.1 GeV" 1 , does only lead to irrelevant 
changes). These are F n = 100 MeV , D = 0.75 , F = 0.5 , C = 1.5 , D R = 0.6 , F R = 0.11 
and we set A = 1 GeV. The various mesonic LECs Ll(X) are the central values taken 
from the compilation of Bijnens et al. in Ref . |j47] . For these central values, all baryon 
masses are fitted exactly, together with cx^a^O) = 45 MeV. The theoretical uncertainties 
induced by the spread of these parameters, in particular due to the ±10 MeV uncertainty 
in ct^jv(O) p2]| , will be discussed in the next subsection. We find for the octet baryon mass 



in the chiral limit using Eq . (|65|) 

m= 767 MeV . (79) 

The quark mass expansion of the baryon masses, in the notation of Eq. (]T6|) , reads 

m N = m (1 + 0.34- 0.35 + 0.24) , 

m A = m (1 + 0.69 - 0.77 + 0.54) , 

ms = m (1 + 0.81 - 0.70 + 0.44) , 

m E = m (l + l.io - 1.16 + 0.78) . (80) 

We observe that there are large cancellations between the second order and the leading 
non-analytic terms of order g 3 , a well-known effect. The fourth order contribution to the 
nucleon mass is fairly small, whereas it is sizeable for the A, the S and the H. This is partly 
due to the small value of m, e.g. for the S the leading term in the quark mass expansion 
gives only about 55% of the physical mass and the second and third order terms cancel 
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almost completely. The fourth order contributions are indeed dominated by the one-loop 
graphs with insertions from as conjectured by Jenkins and Manohar ^ [|. However, 

one can not neglect the terms with insertions from the remaining dimension two terms, 
which are proportional to the 6j and stem from relativistic 1/m corrections. In contrast, 
the contributions from the local terms C^ B are fairly small, i.e. one does not need to know 
the LECs di very accurately. From the chiral expansions exhibited in Eq.(|80"D one can 
not yet draw a final conclusion about the rate of covergence in the three-flavor sector of 
baryon CHPT. Certainly, the breakdown of CHPT claimed in ref.0] is not observed. On 
the other hand, the conjecture [^l|] that only the leading non-analytic corrections (LNAC) 



~ m?J 2 are large and that further terms like the ones ~ m 2 are moderately small, of the 
order of 100 MeV, is not supported by our findings. 

We now turn to the cr-terms and the strangeness content of the nucleon. The pion- 
nucleon cr-term is used in the fitting procedure. It is, however, instructive to disentangle 
the various contributions to 0^(0) of order q 2 , q 3 and g 4 , respectively, 

a wN (0) = 58.3 (1 - 0.56 + 0.33) MeV = 45 MeV , (81) 

which shows a moderate convergence, i.e. the terms of increasing order become succes- 
sively smaller. Still, the q A contribution is important. Also, at this order no tttt rescattering 
effects are included. Rewriting the a-term as §| 

^tv(O) = ^— (82) 

i - y 

we find for the strangeness fraction y and for a 

y = 0.21 , a = 36 MeV . (83) 



The value for y is within the band deduced in ref. ||22|| , y = 0.15±0.10 and the value for 
<t compares favourably with Gasser's estimate, a = 33 ± 5 MeV 0. 

Finally, we consider the kaon-nucleon a-terms and the various scalar form factors. As 
stressed in section |3.3| , there appear undetermined renormalization constants at order g 4 
as long as one works in the isospin limit m u = ma- These are expected to be of order one. 
Indeed, for the pion-nucleon a-term one can calculate this constant (called a n ) since the 
full renormalization has been performed {(3r = 0), 

o-^(O) = (50.1 - 14.4 a w ) MeV -> a n = 0.36 , (84) 

i.e. a n has indeed the expected size. For the kaon-nucleon a-terms, we can give the 
results up to two constants (called axi and a^), 

a$]f(0) = (369 - 306 a K1 ) MeV , 

^Sv(0) = (934 — 437 a^ 2 ) MeV , (85) 
with the respective chiral expansions 

0-^(0) = (528 - 524 + 365 - 306 a K1 ) MeV , 

o-^(O) = (290 - 49 + 693 - 437 a K2 ) MeV , (86) 
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where the terms refer to the orders q 2 , q 3 , q 4 and g 4 , respectively (the g 4 contribution 
independent of the renormalization constant axi is shown separately). These numbers 
agree with the q 3 calculation of ref.Q (to that order). Varying clki and ax2 between 0.5 
and 1, we have 

(7^(0) = 73... 216 MeV , 0^(0) = 493 ... 703 MeV . (87) 

These numbers should be considered indicative and can only be sharpened in a calculation 
with m u ^ TUd- For a discussion of the various a-term shifts to the pertinent Cheng- 
Dashen points, we refer to ref. [30]. 



5.2 Theoretical uncertainties 

In the previous paragraph, we gave the results for the central values of the input param- 
eters. Here, we will discuss the spread of the results due to the uncertainties related to 
these numbers. 

Consider first the dependence on the coupling constant C and the values of the axial- 
vector couplings F and D. For comparison with our central values, we also use \C\ = 1.2 
determined by Butler et al. MM, \C\ = 1.8 from the decay A — ► Nj and D = 0.85 ± 0.06, 



F = 0.52 ± 0.04 given by Luty and White fl49fl . The results depend only very weakly 



on these parameters, i.e. they vary within a few percent. For the case of C, this weak 
dependence stems from the fact that C only changes the value of e& whereas the changes 
in the much more important /3 a and 8& are absorbed in the new fit values of a^ 2 ' 4 '. The 
weak dependence on the actual values of F and D is due to compensating contributions of 
third and fourth order and the already mentioned dominance of the tadpole graphs with 
the symmetry breakers ~ &o,zxf- Consider now variations in the renormalization scale A. 
The latter dependence is induced since we estimate the LECs from resonance exchange 
and would disappear once all LECs could be determined from data. In table 3 we show 
results for the range 0.8 GeV < A < 1.2 GeV, for the central values of F w , F,D,Fr 
and Dr. The strangeness fraction y is most notably affected. The chiral series for the 
masses converges quicker for lower values of A. In table 4, we vary the value of ova^O) 
generously by ±10 MeV. Again, the strangeness fraction shows the largest variation. All 
these variations are essential linear in Scr n ]^(0). Finally, we remark that varying the 
pseudeoscalar decay constant F w between 93 and 113 MeV does also not alter any of the 
previous numbers drastically. We therefore assign the following theoretical uncertainty to 
the results for the average octet mass in the chiral limit, the strangeness fraction y and 
a, in order 

m= 767 ± 110 MeV , y = 0.21 ± 0.20 , a = 36 ± 7 MeV . (88) 

These uncertainties do not include the possible effects of higher orders which can only 
be assesed if one performs such a calculation. As an indication of genuine two loop 
contributions, we quote the results of the diagrams evaluated in paragraph [R| We find 
(setting A = 1 GeV) 

m % a) = -52 MeV m% f) = -13 MeV , (89) 
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which are individually very small, i.e. well within the uncertainties discussed above. It is 
expected that the major contribution at order q 5 does indeed come from one-loop graphs 

(3) 

with insertions from and not from the genuine two loop diagrams. To quantify this 
statement, such an order q 5 calculation has to be performed. That, however, goes beyond 
the scope of the present paper. 

6 Summary and outlook 

In this paper, we have considered the chiral expansion of the groundstate baryon octet 
baryon masses and the pion-nucleon cr-term to quadratic (fourth) order in the quark 
(Goldstone boson) masses, in the framework of heavy baryon chiral perturbation theory. 
The pertinent results of this investiagtion can be summarized as follows: 

o We have constructed the most general effective Lagrangian to fourth order in the 
small parameter q (external momentum or meson mass) necessary to investigate 
the scalar sector. Besides the standard dimension two symmetry-breaking terms, 
Eq.(||), it contains further dimension two operators with derivatives acting on the 
Goldstone boson fields and kinematical 1/m corrections, cf. Eq.(|9]) and Eq.(|T3|). 

o We have given the complete expressions for the baryon masses (^at,a,s,h) at order 
q A togther with the pion-nucleon cr-term. At this order, divergences appear (in 
contrast to 0(q 3 ), where the one loop corrections to the masses and ova^O) are 
finite). The renormalization procedure to render the baryon self-energies and also 
cVa^O) finite (by use of the Feynman-Hellmanm theorem) involves additional terms 
as listed in Eq. (p8|) . The renormalization of the kaon-nucleon cr-terms and the 
corresponding scalar form factors is further complicated by additional momentum- 
independent divergences, as detailed in appendix |B|. 

o There appear seven low-energy constants (LECs) at order q 2 (called hi) and seven 
more at order q 4 (called di) to calculate the masses and pion-nucleon cr-term. Two 
of the latter amount to quark mass renormalizations of two of the hi. Since there 
do not exist enough data to fix all these, we have estimated them by means of 
resonance exchange. Besides standard tree graphs with scalar meson exchange, 
this involves Goldstone boson loops with intermediate baryon resonances (spin- 
3/2 decuplet and the spin-1/2 (Roper) octet) for the scalar-isoscalar LECs. We 
have discussed a consistent scheme how to implement resonance exchange under 
such circumstances, i.e. which avoids double-counting and abids to the strictures 
from analyticity. Within the one-loop approximation and to leading order in the 
resonance masses, the analytic pieces of the pertinent graphs are still divergent, 
i.e. one is left with three a priori undetermined renormalization constants (/?a, 
and (3r). These have to be determined togther with the finite scalar couplings F s 
and Ds and the octet mass in the chiral limit. Using the baryon masses and the 
value of ova^O) as input, we can determine all LECs in terms of one parameter, (3r. 
We derive a bound on this parameter and show that observables are insensitive to 
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variotions of (3r within its allowed range. We have also demonstrated that the effects 
of two (and higher) loop diagrams can almost entirely be absorbed in a redinition 
of the one loop renormalization parameters. 

o Within this scheme of estimating the LECs we determine the baryon mass in the 
chiral limit, denoted m, m= 770 ± 110 MeV (accounting also for the uncertainty in 
certain input parameters like e.g. the pion-nucleon cr-term). For the strangeness 
fraction y we find y = 0.21 ± 0.20, consistent with dispersion-theoretical determi- 
nations. This translates into a = 36 ± 7 MeV, compare Eq . fl82|) , in good agreement 
with a previous calculation 0]. 

o The chiral expansions for the nucleon mass and the pion-nucleon a-term are moder- 
ately well behaved whereas for the hyperon masses m^j; = there still appear sizeable 
corrections at order q 4 . This is partly due to the almost complete cancelations be- 
tween the terms of order m q and m^J 2 and the smallness of the baryon mass in 
the chiral limit. A definite statement about the convergence of three-flavor baryon 
CHPT can thus not yet be made. 

o We have also estimated the two kaon-nucleon a-terms, which take the ranges given 
in Eq. (|86|) based on dimensional analysis for the appearing renormalization constants 
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A Baryon mass coefficients at order q A 

Here, we give the coefficients with i — 1, 2, 3, 4 and P,Q = n, K,rj for the A, the £ 
and the S. All di are understood as <i[(A). 
A : 



Hi a = -4(44 + ^ 4 + d 7 + 3d 8 ) , ef A = -16(-d 3 + ^d 4 + d 7 + d 8 ) 



2 

el K h = 16 (- rf 7 + ^8 + |d 3 + U±) , e^ A = 2b D + 6b - 46 3 - 66 8 - ^- 

20 4 D 2 -I- 9 F 2 

4a = irb D + 8b - 126a - -6 3 - 86 8 



3 3 3 m 

D 2 

3 m 



: 3jA 26 - 46 3 - 26 8 , eV A = -16 D 2 b D , e\\ = 16 D 2 b D 
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<a = -(D 2 + 9F 2 )(^b d -2b )-(D-3F) 2 (-2b F + b ) 

-(D + 3F) 2 (2b F + b ) , 



16, 

-2(D + 3F) 2 (b + b D - b F ) , 



<a = P 2 + 9F 2 )^6 d + 46o)-2( j D-3F) 2 (6 + 6d + M 



= -4(44 + +d 7 + 3rf 8 ) , eg = 16(-d 7 + 4) , ef E = "16(4 + d 8 ) , 

D 2 + 6F 2 

e - E = 66 D + 66 - 8&i - 46 3 - 66 8 5 , 

3 m 

6 m 

4 D 2 Ifi 

eJ E = 2& --6 3 -2& 8 -— ,eg = ^D 2 & D , 
<J 3 m <J 

<e = (^ 2 + 4F 2 )(66 D + 36 )-4 J D 2 (-^z ? + ^o)-12(26 I) + 6 ) , 

eg = Q( D 2 + F 2)(2b D + b )-3(D + F) 2 (-2b F + b )-3(D-F) 2 (b + 2b f ) , 
ef« = 12(D 2 + F 2 )b -6(D + F) 2 (b + b D + b F )-6(D-F) 2 (b + b D -b F ) , 

<J = Ibn ,e:% = . (A.2) 

ej >s = -4(4di - 2d 5 + d 7 + 3d 8 ) , efjg = -16(di + d2 + d 3 + d 5 + d 7 + d 8 ) , 
eg = 8(4di + 2d 2 + d 5 - 2d 7 + 2d 8 ) , 

£3,3 = 3(6 D - b F + 2b -b 1 + b 2 -h- 2b 8 ) - 3 (D ~ , 

4 m 

* . , . „, , „, , 5D 2 + 6DF + 9F 2 



4~ = 2(3b D + b F + Ab - 3h - b 2 - 3b 3 - Ab 8 ) - 



6 m 



e|j s = 2b -3b 1 -b 2 -h 3 -2b 8 - {D + 3 p 2 ,ell = , e\l = , 

<J 12 m 

<s = ( 5 ^ 2 + + 9i?2 ) (- 26 f + 60) " I (D + F) 2 (2b D + b ) 

- l -{D-3F) 2 {3bo-2b D ) , 
eff = 2(5D 2 + 6DF + 9F 2 )(b D + b F + b )-9(D + F) 2 b 

~(D-3F) 2 (3b + 4b D ) , 

= -b D -\bp ,efa = \{b D + b F ) . (A.3) 
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B Renormalization of the kaon— nucleon a— terms and 
scalar form factors 



Here, we discuss the renormalization related to the two kaon-nucleon cr-terms and the 
three scalar form factors a n N(t) and a K ' N (t). Since we are dealing with composite op- 
erators, we use the standard procedure and generalize the original Lagrangian used to 
calculate the baryon masses, C = X^i A> m the following way 



, dm q 



£ -> £ + J2' c Q,i ~l ' Q = u,d,s } (B.l) 



where the c q j are sources and the prime on the sum indicates that not all values of q and 
i are necessarily taken. Consider now a diagram labelled 'j', which leads to a divergent 
contribution to the cr-term. The renormalization for that graph proceeds via 

dC, dCj 1 , ^/ „ / dC 



^LY' f, k J-^±-) (B.2) 



dm q dm q ^ kql v „„ Uq . 

with the f q \k,j are Clebsch-Gordan coefficients. To order g 4 , there are two types of 
divergences, namely ~ LM 2 t and ~ LM 4 . The t-dependent ones can be renormalized 
with the help of two counter terms of the form 

A = Tr(B{ X+ , [£>„, [D",B]]}) , A = Tr(B[ X+ , [D^B]))) . (B.3) 

For the t-independent ones, one needs six independent terms, 

A = Tr(B[ X+ ,[ X+ ,B}\) , C A = Tr(B{ X+ Ax + ,B]}) , £ 5 = Tr(B{ X+ ,{ X+ ,B}}) 

A = Tr(5 X+ ) Tr( X+ B) , A = Tr(B\x+, B)) Tr( X+ ) , A = Tr(BB) Tr( X+ ) Tr( X+ ) . 

(B.4) 

Consider as an example the composite operator uu and the corresponding matrix element 
< p'\uu\p >. A specific contribution to it is given by the diagram in fig. 2 with an insertion 
of the &i-vertex, leading to (cf. Eq.(|31|)), and neglecting for simplicity the n° — rj mixing, 

I a (t) = i (AMI -t)L + finite terms (B.5) 

Expanding the operator uu and the &i-term of the effective Lagrangian to second order 
in the meson fields, the divergent parts from l a contribute as follows 

< p'\uu\p > div = B L (*-Ml + 2M 2 K + ±AdA + ^LBLt. (B.6) 

7T 7T 

The t-dependent divergence is renormalized with the help of C\ via 

2b, r S£, ac, aCn 

-Jf L {- 3 a^- 2 d^-d^J (B7) 
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i.e . for the Clebsch-Gordan coefficients 

fu,i,j = -6 61 , = -4 61 , = -6 b ± , (B.8) 

and zero otherwise. The t-independent divergences in the matrix element < p'\uu\p > 
are renormalized with the help of £3, £5, Cq and 



I 12<9m u 12dm d 8 dm u 8 dm d 8dm s dm u 2dm u 
with the Clebsch-Gordan coefficients 

/«,3,j = — 24^ x ' •^ d ' 3 'j' = ~ 24^ ' ^ U ' 5 ' J, = ~~ 16 ^ ' ^' rf ' 5 ' : ' = ~~ 16 ^ ' 

1 1 
fs,5,j = ~~[Q bl ' ^ U,6,j = ^ ' ^ u ' 8,j = ~4 &1 ' (B.10) 

C (J-kN^) and strangeness fraction coefficients at 0(q A ) 

Here, we collect the coefficients and efy® with i = 1, 2, 3 and P,Q — ir, K, r\. 
Pion-nucleon a-term: 

TrK 1 7T , 7T7T , ^ 7T_ft" , ^ WT 

2,7V e 3,7V e 4,7V "r 2 4,7V g~ e 4,7V 5 
if , 1 | wK | 1 «tt 





— ze l,7V + 2 




— e l,7V + e 2,7V 




2 „ 1 
= 3 e ^ + 3 f 


€ 2,a 


~~ fc 3,7V ^ fc 4,7V 


K-k 

e 3,cr 


— e 4,7V > e 3,cr 



2 3,7V 2 4 ,7V 4 ,7V <j C 4,jV 5 

_|_ I^W ^ — 9^ _|_ Of 71 " 71 " _|_ I, 
+ 3 4,/V ' e 2,cr — ze 3,7V + ze 4,7V + 2 



9 9 1 

KK v — C V i ,W i ,>pr i 
4,7V ) e 2,<r — 3 e 3,7V g~ e 4,7V e 4,7V ^ e 4,7V 5 

1 JOr ^ _ 1 rpr Kr; _ 1 r]K /p -, \ 

2 4,7V 5 e 3,<r — 3 e 4,7V ) e 3,o- — 3 4,7V ■ V^- 1 ) 

Notice that the 6(2,3)0- and the e( 3i4 )7v include the factor 1/A£ which is not made explicit 
here. 

Strangeness fraction: 



(C.2) 



f 77 — 


1 kK 


+ 


1 Ktv 

2 4 > N 


2 

+ 3 4,7V 


> e l,j/ - e l,7V + 


-e K + 

2 e 3,7V + 


2 „a- 
2 e 4,/v ? 




e hy - 


3 €l ' N 


+ 


3 e 3 ,7V 




U > e 2,y — e 3,7V 








V 

e 2,y - 


3 3 > N 


+ 


1 nK 
2 e 4,7V 


c -kK 
' C 3,j/ - 


2 e 4,7V > e 3,y ~ 


U ' € 3,y 


9 

— 3 4,/V ) 


Kr, _ % V K 
e 3,j/ - 3 e 4,7V 
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D Determination of the Roper— octet coupling con- 
stants Dr and Fr 

The only listed deacys to determine the coupling constants Dr and Fr are the A/"* (1440) — > 
Ntc and the A*(1600) — > £71". Consider first the A* decay. The width follows via 

wj\t\ 2 , (d.i) 



57rm 



with | q % \ = 334 MeV the pion momentum in the restframe of the A* and we use the 
spinor normalization u{p)u{p) = 2mR. Straightforward calculation gives 



2E V (E s E n + ql) - Ml (Ez + m s ) 



167rm\*F^ 

= 149 MeV D 2 R , (D.2) 

with E n = (Ml + ql) 1 ' 2 and E s = (m| + g^) 1/2 . From the PDG value T(A* -> Stt) = 
(150 ± 100) MeV • (0.35 ± 0.25) = 52.5 MeV we derive 

D R = 0.60 ± 0.41 . (D.3) 

The sign of Dr is, of course, not fixed. We will chose it to be positive in accordance with 
the ground state octet D coupling. The Af*(1440)A% effective Lagrangian is 

C>n*n-k = VR^ml^u^^N + h.c. , 



g A VR = Dr + Fr . (D.4) 

Here, \/R = 0.53 ± 0.04 and g& = 1.33 as given by the Goldberger-Treiman relation. As 
explained in some detail in Ref . [fll]] , we use the width of the N* determined from the speed 
plot, not the the model-depedent Breit-Wigner fits, r to t = 160 ± 40 MeV. Assuming Fr 
to be positive (as is the equivalent hyperon coupling), we find by using Eq.( 



^ = 0.11^=0.41 . (D.5) 
These are the values of the Roper couplings used in the main text. 
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Tables 



B 


Xb [GeV 2 ] 


D B [GeV 2 ] 


D 5 B [GeV 4 ] 


D% [GeV 2 ] 


N 


0.030 


(-1.165,+0.233) 


(-0.123,-0.079) 


(0.114,0.029) 


A 


0.065 


(-2.013,-0.783) 


(+0.077,-0.053) 


(0.219,0.035) 


£ 


0.058 


(-3.626,-0.631) 


(-0.802,-0.186) 


(0.506,0.058) 




0.096 


(-3.205,-1.648) 


(-0.390,-0.108) 


(0.471,0.052) 



Table 1: Numerical values of the state-dependent coefficients in Eq.(^). The first and 
second term in the brackets refers to the decuplet and Roper-octet contribution, respec- 
tively. 





b 


b D 


b F 


h 


b 2 


h 


h 


(3r = 


-0.606 


0.079 


-0.316 


-0.004 


-0.187 


+0.018 


-0.109 


Ref.[42] 


-0.493 


0.066 


-0.213 


+0.044 


-0.145 


-0.054 


-0.165 


Ref.0 


-0.750 


0.016 


-0.553 











Table 2: Low-energy constants from C\ B in GeV 1 for (3 R = 0. 
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A 


7 


7 


7 

b F 


O 

m 


(2) 

my 


(4) 

my 


a 


y 


[GeV] 


[GeV L \ 


[GeV X J 


[GeV X J 


[MeVJ 


[MeVJ 


[~\ T T T~\ 

[MeVJ 


[MeVJ 




0.8 


-0.643 


0.085 


-0.369 


849. 


244. 


119. 


49.4 


-0.10 


1.0 


-0.606 


0.079 


-0.316 


767. 


261. 


183. 


35.6 


+0.21 


1.2 


-0.574 


0.072 


-0.281 


613. 


266. 


232. 


25.9 


+0.42 



Table 3: Theoretical uncertainties due to the renormalization scale A. For comparison, 
we give the second and fourth order contribution to the nucleon mass, m$ and m$ , 
respectively. The third order contribution is m$ = — 272 MeV. We set (3 R = 0. 



<t*n(P) 


b 


b D 


b F 


o 

m 


(2) 

my 


(4) 

my 


a 


y 


[MeV] 


[GeV- 1 ] 


[GeV- 1 ] 


[GeV- 1 ] 


[MeV] 


[MeV] 


[MeV] 


[MeV] 




35 


-0.508 


0.082 


-0.319 


910. 


163. 


139. 


37.4 


-0.07 


45 


-0.606 


0.079 


-0.316 


767. 


261. 


183. 


35.6 


+0.21 


55 


-0.691 


0.074 


-0.310 


625. 


360. 


226. 


37.1 


+0.34 



Table 4: Theoretical uncertainties due to the pion-nucleon cr-term. For comparison, 
we give the second and fourth order contribution to the nucleon mass, m$ and mff , 
respectively. The third order contribution is ruff = —272 MeV. We set (3r = 0. 



Figure captions 

Fig.l One-loop graphs with exactly one insertion from C^ B (circlecross) . Goldstone boson 
renormalizations are not shown. 

Fig. 2 Contribution to the scalar form factor. The double-line denotes the insertion of 
m{uu + dd). 

Fig. 3 Baryon resonance excitation involving pion loops. The double-line represents the 
decuplet or the even-parity Roper octet. Solid and dashed lines denote the ground 
state octet baryons and Goldstone boson fields, respectively. 
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Fig.4 Scalar meson excitation. The double-line represents the scalars and the circlecross 
a quark mass insertion ~ x+- 

Fig. 5 Baryon resonance excitation: Two-loop graphs. The double-line represents the 
decuplet or the even-parity Roper octet or a combination thereof. Solid and dashed 
lines denote the ground state octet baryons and Goldstone boson fields, respectively. 
Graphs of the same topologies with one or two groundstate propagators are not 
shown. 
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Figure 5 



41 



